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Abstract. We discuss our numerical results on the properties of the 5 = 1/2 frustrated Ji- 
Ji Heisenberg model on a square lattice as a function of temperature and frustration angle 
4> = tan - (J2/J1) in an applied magnetic field. We cover the full phase diagram of the model 
in the range — ir < (j> < n. The discussion includes the parameter dependence of the saturation 
field itself, and addresses the instabilities associated with it. We also discuss the magnetocaloric 
effect of the model and show how it can be used to uniquely determine the effective interaction 
constants of the compounds which were investigated experimentally. 



1. Introduction 

The ground-state and low-temperature properties of low-dimensional frustrated spin systems 
are of great interest due to their unusual and sometimes unexpected temperature and magnetic- 
field dependence. The origin for this lies in the enhanced quantum fluctuations which determine 
their behaviour. A prime example for such a system is the spin-1/2 frustrated J\-Ji Heisenberg 
model on a square lattice, being one of the most intensively studied models for spin systems. It 
has a Hamiltonian of the form 

H = J 1 J2 SiSj- + J 2 £ SiSj, (1) 

(y>i (*i> 2 

where J\ denotes the exchange interaction between neighbouring spins, and J 2 labels the 
exchange interaction along the diagonals between next nearest neighbours. 

In order to facilitate the description of the model and its thermodynamics, it is convenient to 
introduce two parameters to characterise it, namely an overall energy scale J c and a frustration 
angle 4> defined by 



J c = ^J2 + J| an d cj) = tan -1 (J 2 /Ji) , (2) 

such that Ji = J c cos 4>, and J 2 = J c sin </>. Except where stated otherwise, the numerical results 
discussed in the following are all obtained using the finite-temperature Lanczos method [I; on a 
20-site square with periodic boundary conditions. 



2. Phase diagram and experiments 

Figure n shows the phase diagram of the model, together with some extra lines referring to the 
compounds which were investigated experimentally. At a classical level, the phase diagram has 




J2/J1= -oo 



Figure 1. Schematic phase diagram of the J\-Ji model. The boundary between FM and NAF 
phase is the line J\ = 0, J2 < 0. Numbers are ratios of exchange couplings J2/J1 where zero 
point fluctuations destroy the relevant order parameter ^ The dashed lines labeled by 
the chemical formulae for the experimentally known "J1-J2 compounds" are lines of constant 
@CW = (Ji + ^2)Ab as determined from H = heat capacity and magnetic susceptibility 
measurements 011]. 

three main phases: For —tt/2 < <f) < tan _1 (l/2), the Neel antiferromagnet (NAF) characterised 
by an ordering vector Q = (vr,7r), a ferromagnetic phase (FM) for tt — tan~ 1 (l/2) < <fi < —it/2, 
and a collinear phase (CAF) with Q = (ir,Q) or (0,7r) for tan _1 (l/2) < < tt - tan -1 (1/2). 
At the boundary between CAF and NAF, the Ising domain wall energy between the states 
vanishes, and arbitrary stripes of the two phases can be formed. At the CAF/FM boundary, 
coplanar spiral states with q = (2im/m, 0), where {n,m} are integers, all become degenerate. 
Quantum fluctuations lift these degeneracies and lead to disordered regions of finite width at 
these two boundaries having continuous classical degeneracies. A more detailed discussion of 
the properties of the model at zero magnetic field is given in References |2j and the references 
cited therein. 

The J\-Ji model is appropriate for the description of a number of quasi two-dimensional 
vanadium compounds. They can be grouped in two families of type Li2VOX04 (X = Si, Ge), 
and AA'VO(P0 4 ) 2 (A, A' = Pb, Zn, Sr, Ba), respectively [SHU. A11 compounds feature V 4+ 
ions with S = 1/2 surrounded by oxygen polyhedra, forming layers of J1-J2 square lattices 
which are weakly coupled in the third spatial dimension. 

Thermodynamic coefficients have been measured to determine the effective interaction 
constants jSH3]. For zero magnetic field, this is not completely possible. An ambiguity remains 
when fitting the experimental xC^O an d Cy(T) using a high-temperature series expansion: We 
always have two possible points <fi = <f>± in the phase diagram symmetric around <fi = ir/4. This 
situation does not change qualitatively when including higher-order terms in the expansion used 
for the fits [3]- 

3. Finite fields 

Fortunately, all J\-J% compounds known so far have values J c = (D(10K), which corresponds to 
a magnetic-field range accessible in pulsed high-magnetic field experiments. 

We have calculated the saturation field of the magnetisation as a function of the frustration 
angle 4> and the magnetic field H. The left-hand side of Figure [2] shows the result for a 20-site 




Figure 2. Left: One-magnon instability (solid line) and two-magnon instability (dashed line) 
field H of the fully polarized state of the J1-J2 model on a 20-site square as a function of 
the frustration angle <p = tan" 1 (J2/ J\) ■ The one-magnon instability field corresponds to the 
saturation field of the magnetisation obtained by linear spinwave theory. The energy unit is 

J c = \J J I + J 2 ■ Right: Difference between the AS = 1 and AS = 2 instabilities as a function 
of the frustration angle for a cluster with 16 sites (dotted line), 20 sites (solid line), and 24 sites 
(open circles). 



cluster with periodic boundary conditions as well as the saturation field obtained from linear 
spinwave theory. We define the k th saturation field as the instability of the fully polarised 
(saturated) state against an excitation with k magnons, 
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where E denotes the ground-state energy per site of the cluster of size N, and S^ ot is the 
component of the total spin of the system in the direction of the magnetic field H. In 
this notation, the thermodynamic saturation field of the magnetisation is given by -ff S at 
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The result from linear spinwave theory for the infinitely large system reads 
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with z = 4 and 5 = 1/2, and is displayed as solid line in FigureEl(left). It is in perfect agreement 
with the numerical result for hQ. 

For values —n/2 < eft < tt/2 [J\ > 0), the saturation field is determined by the A5 = 1 
instability. Inside the collinear phase, at the point eft = tt/2, J\ vanishes, and the lattice 
decouples into two sublattices. Therefore, the AS = 1 and AS = 2 spin flip energies must 
have the same value. For ferromagnetic J±, for the finite-size clusters considered, the AS = 2 
instability determines the saturation field until the point <j> ~ 0.97T where the ferromagnetic phase 
is reached. This is seen from Figure 12 (right), where we have plotted the difference H^ t — H^ t 



for three different cluster sizes. While the overall size of H^ t — H^l decreases monotonically 
with the size of the cluster, at least in the NAF and CAF phases, the crossover described above 
happens always exactly at J\ = independent of the cluster size. 
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4. Magnetocaloric effect 

The magnetocaloric effect (MCE) is the adiabatic temperature change of a sample by applying 
a magnetic field, given by 
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Figure 3. Left: H-cj> contour plot of the MCE of the J\- J2 model at T = 0.2J c //cb on a 20-site 
square with periodic boundary conditions. The calculation was done on a grid of 330 x 200 points. 
Right: Field dependence of the MCE for those two particular values of eft = <p± determined for the 
compound SrZnVO(P04)2 [2], again at a temperature T = 0.2J c /£;b- Note the large difference 
of the maximum positions. 



f8T\ ^_(0S_\ /fdS\ T{H) / dM(H, T) \ 

\dHJs \dHj T /\dTj H ~ C V (H,T)\ dT ) H > {> 

where S is the entropy, and M the magnetisation. We have calculated the MCE on a 20-site 
square as a function of T, H, and (j). Figure |3] (left) shows the MCE for fixed temperature 
T = Q.2J c /k-Q. Its main feature is the pronounced anomaly near the saturation field where 
the fully polarised state is reached, corresponding to the sudden entropy change taking place 
upon reaching saturation. The saturation magnetisation is reached slightly above the maximum 
position of the MCE. The oscillating structure of the MCE below -ff S at m the CAF phase is a 
finite-size effect which we expect to vanish in the thermodynamic limit. Note that in the Neel as 
well as in the collinear phase, the MCE is for some range of fields negative, indicating that the 
sample will cool down with increasing field, opposite to what is expected for simple paramagnets. 

The enhanced quantum fluctuations due to the frustration are evidenced by an enhanced 
magnetocaloric effect in the disordered region between the CAF and the FM phases, and close 
to the region between the CAF and the NAF phases. 

The right-hand side of Figure El shows two vertical cuts through the left-hand side plot for 
those values = <f>± which have been determined experimentally from zero-field heat capacity 
and susceptibility measurements of the compound SrZnVO(P04)2 3 . In contrast to the zero- 
field data, the two maximum positions of the MCE anomalies differ by roughly a factor two. 

The maximum saturation field of the J1-J2 model is reached at ^ max /7r = tan _1 (2)/7r ~ 0.35 
with gfiftH™^ = 4.5J C . Since the available compounds all have values J c = O(10 K), their 
saturation fields have values H sllt < 40 T, assuming g = 2. These fields are accessible; therefore, 
a measurement of the field dependence of the magnetocaloric effect can help to clarify the nature 
of the ground state of the J\-Ji compounds discussed here. 
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